On stable determination of potential by boundary measurements 



R. G. Novikov 



CNRS, Laboratoire de Mathematiques Jean Leray (UMR 6629), Universite de Nantes, 
BP 92208, F-44322, Nantes cedex 03, France 
t*"^ ' e-mail: novikov@math.univ-nantes.fr 

o" 

Abstract. We give new stability estimates for the Gel' fand-C alder on inverse bound- 
ary value problem. 

1. Introduction 

Consider the equation 



PL, 



-A^ + v(x)i/) = 0, xeD, (1.1) 



where 



J ■ D is an open bounded domain in 

dDeC 2 , veL°°(D), d>2. 



(1.2) 



We assume also that 



is not a Dirichlet eigenvalue for 

(1.3) 

the operator — A + v in D. 



> 

\ Equation (1.1) arises, in particular, in quantum mechanics, acoustics, electrodynamics 

Formally, (1.1) looks as the Schrodinger equation with potential v at zero energy. 
CN ■ Consider the map $ such that 

— I =Q(iI>\ ) (14) 

Q v \dD *W\ dD ) ^ > 

for all sufficiently regular solutions ip of (1.1) in D = D U dD, where v is the outward 
normal to dD. The map $ is called the Dirichlet-to-Neumann map for equation (1.1) and 
is considered as boundary measurements for (1.1). 

We consider the following inverse boundary value problem for equation (1.1): 

Problem 1.1. Given $, find v. 

This problem can be considered as the Gel'fand inverse boundary value problem for the 
Schrodinger equation at zero energy (see [G], [Nol]). This problem can be also considered 
as a generalization of the Calderon problem of the electrical impedance tomography (see 
[C], [SU], [Nol]). Concerning results given in the literature on Problem 1.1 (in its Calderon 
or Gel'fand form ) see [KV], [SU], [HN] (note added in proof), [Nol], [Al], [Nal], [Na2], 
[BU], [P], [Ma], [No2], [No4], [Am] and references therein. 

In the present article we show that the Alessandrini stability estimates of [Al] for 
Problem 1.1 in dimension d > 3 (see Theorem 2.1 of the next section) admit some princi- 
ple improvement. Our new stability estimates (see Theorem 2.2 of the next section) are 
obtained by methods developed in [No2], [No3], [No4]. These methods include, in partic- 
ular: (1) the d- approach to inverse "scattering" at zero energy in dimension d > 3, going 



R.G. Novikov 



back to [BC], [HN], and (2) the reduction of Problem 1.1 to inverse "scattering" at zero 
energy, going back to [Nol]. 

The present article is organized as follows. In Section 2, we formulate and discuss old 
and new stability estimate for Problem 1.1. In Section 3, we remind (a) definition and 
some properties of the Faddeev functions and (b) formulation of the inverse "scattering" 
problem for the Schrodinger equation at zero energy (Problem 3.1). In Section 4, we 
remind formulas and equations of [Nol], [No2] reducing Problem 1.1 to Problem 3.1. In 
Section 5, we remind an approximate reconstruction method of [Nol] for Problem 1.1. In 
Section 6 we prove Theorem 2.2 in the Born approximation. 

2. Stability estimates 

We assume for simplicity that 



D is an open bounded domain in R d , dD e C 2 , 

v e W m,1 (R d ) for some m > d, suppv C D, d>2, 



(2.1) 



where 
where 

Let 
Let 



W m,1 (R d ) = {v : d J v e L 1 (M d ), |J| < m}, mGNUO, (2.2) 
Je(NU0) rf , |J| = X>, »>c,r) J," r( ;',' 



i=i 



dxi 1 . . . dx J d a 



i= max \\d J v\\ Tl ,-rad.. (2.3) 

\J\<m L ( M ) 



(2.4) 



|| A || denote the norm of an operator 
A : L°°(dD) -> L°°(dD). 

We remind that if v±, V2 are potentials satisfying (1.2), (1.3), where D is fixed, then 

$i — $2 is a compact operator in L°°(dD), (2.5) 

where $i, $2 are the DtN maps for v±, V2 respectively, see [Nol], [No2]. Note also that 
(2.1) =► (1.2). 

Theorem 2.1 (variation of the result of [Al]). Let conditions (1.3), (2.1) hold for 
potentials v\ and v^, where D is fixed, d > 3. Let \\vj\\ m! i < R, j = 1, 2, for some R > 0. 
Let $1, $2 denote the DtN maps for v\, v%, respectively. Then 

\\vi - v 2 \\ L °° (D ) < Ci(ln(l + ||$! - ^ir 1 )) -01 , (2.6) 

where C\ = C\(R, D,m), ol\ = (m — d)/m, ||$i — $2!) is defined according to (2.4). 
Theorem 2.1 follows from formulas (3.9)-(3.11), (4.1) (of Sections 3 and 4). 
A disadvantage of estimate (2.6) is that 

cti < 1 for any m > d even if m is very great. (2.7) 
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Theorem 2.2. Let the assumptions of Theorem 2.1 hold. Then 

\\vi - v 2 \\l°°{d) < C 2 (ln(l + ||$! - ^ir 1 ))— 2 , (2.8) 

where C 2 = C 2 (R, D,m), a 2 = m — d, ||$i — <& 2 \\ is defined according to (2.4). 
A principal advantage of estimate (2.8) in comparison with (2.6) is that 

a 2 — > +oo as m — > +oo, (2-9) 

in contrast with (2.7). 

In the Born approximation, that is in the linear approximation near zero potential, 
Theorem 2.2 is proved in Section 6. 

For sufficiently small R in dimension d = 3, Theorem 2.2 follows from (3.9) (of Section 
3) and results of [No2] , [No4] . The scheme of our proof for this case is, actually, similar to 
the scheme of our proof for the case of the Born approximation. The main difference is that 
instead of the inverse Fourier transform (used in Section 6) we use now the zero-energy 
inverse " backscattering" transform of [No4] . We plan to give this " nonlinear" " small- norm" 
proof in a separate article. 

In the general case, the proof of Theorem 2.2 is not completed yet. However, except 
restrictions in time, we see no difficulties for completing this proof by methods of [No2] , 
[No3], [No4]. 

3. Faddeev functions 

We consider the Faddeev functions G, ip and h (see [Fl], [F2], [HN], [No2]): 

G{x,k) = e**g{x,k), g{x,k) = -(±) d f Rd 0-1) 

^ k ) = e ikx + f G {x-y, k)v(y)iP(y, k)dy, (3.2) 
JR d 



where x G M d , fceS, 



where (fc, /) G ©, 
We remind that: 



E = {keC d , k 2 = kf + ... + k 2 d = 0}; (3.3) 
Kk, I) = {^Y J Rd e- Ux v(x)^(x, k)dx, (3.4) 

= {k G E, / G E : Imk = Iml}. (3.5) 



AG(x,k) = 8(x), x G M d , fceS; (3.6) 
formula (3.2) at fixed k is considered as an equation for 

t(j = e lkx ^(x,k), (3.7) 
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where \i is sought in L°°(R d ); as a corollary of (3.2), ip satisfies (1.1); h of (3.4) is a 
generalized "scattering" amplitude in the complex domain at zero energy. 

Note that, actually, G, ip, h of (3.1)-(3.5) are zero energy restrictions of functions 
introduced by Faddeev as extentions to the complex domain of some functions of the 
classical scattering theory for the Schrodinger equation at positive energies. In addition, 
G, ip, h in their zero energy restriction were considered for the first time in [BC]. The 
Faddeev functions G, ip, h were, actually, rediscovered in [BC]. 

We remind also that, under the assumptions of Theorem 2.1 

p(x,k) — > 1 as \Imk\ — > oo (uniformly in x) (3.8) 

and, for any c > 1, 

\p{x, k)\ < c for \Im k\ > p±(R, D, m, c), (3.9) 

where x G R d , fc£S; 

v(p) = lim h(k,l) for any peR d , (3.10) 

(k,l)e&, k-l=p 
\Im k\ = \Im l\— >oo 

\v(p) - h(k, I) | < C ^ D ^ R2 f or (fc, Z) G 0, p = fc - /, n 



| Jm k\ = \Iml\ = p > P2(R, D,m), 



where 



v(p) = e ipx v(x)dx, p G R d . (3.12) 



Results of the type (3.8), (3.9) go back to [BC]. Results of the type (3.10), (3.11) (with 
less precise right-hand side in (3.11)) go back to [HN]. Estimates (3.8), (3.11) are related 
also with some important L 2 -estimate going back to [SU] on the Green function g of (3.1). 

For more information on properties of the Faddeev functions G, ip, h, see [HN], [No2], 
[No4] and references therein. 

In the next section we remind that Problem 1.1 (of Introduction) admits a reduction 
to the following inverse "scattering" problem: 

Problem 3.1. Given h on G, find v on R d . 

4. Reduction of [Nol], [No2] 

Let conditions (1.2), (1.3) hold for potentials v\ and v%, where D is fixed. Let ipi, 
hi denote the DtN map $ and the Faddeev functions ip, h for v = Vi, i = 1,2. Let also 
$>i(x, y) denote the Schwartz kernel $>(x, y) of the integral operator $ for v = Vi, i = 1, 2. 
Then (see [No2] for details): 

h 2 (k,i)-h 1 (k,i) = (±-) d J J Mx,-i)(*2-*i)&y)My,k)dydx, (4.1) 

dWD 
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where (k, I) G ©; 

if> 2 (x,k) = if>i(x,k)+ A(x,y,k)^2(y,k)dy, x G 3D, (4.2a) 

J 3D 

A(x,y,k)= R 1 (x,z,k){^ 2 -^i){z,y)dz, x,yedD, (4.26) 

J 3D 

RAx,y,k) = G(x-y,k)+ / G(x - z, k)v 1 (z)R 1 (z, y, k)dz, x, y G M d , (4.3) 

m d 

where feeS. Note that: (4.1) is an explicit formula, (4.2a) is considered as an equation 
for finding ip2 on 3D from ipi on 3D and A on 3D x <9D for each fixed k, (4.2b) is an 
explicit formula, (4.3) is an equation for finding Ri from G and v±, where G is the function 
of (3.1). 

Note that formulas and equations (4.1)-(4.3) for v\ = were given in [Nol] (see also 
[HN] (Note added in proof), [Nal], [Na2]). In this case h x = 0, V>i = e ikx , i?i = AO- 
Formulas and equations (4.1)-(4.3) for the general case were given in [No2]. 

Formulas and equations (4.1)-(4.3) with fixed background potential v\ reduce Problem 
1.1 (of Introduction) to Problem 3.1 (of Section 3). 

5. Reconstruction of [Nol] in the Born approximation 

In the Born approximation, that is in the linear approximation near zero potential, 
we have that 

h(k,l)^v{k-l), (5.1) 

dWD 

where (k, I) e 0, v is defined by (3.12), $o denotes the DtN map for v = 0. 

Formulas (5.1), (5.2) follow from (3.1)-(3.4) and (4.1). Formulas (5.1), (5.2) imply, in 
particular, that 

v(p) « {^) d J J e- il ^ x (^-^ )(x,y)e ik ^ydxdy, (5.3) 

dWD 

k(p) = \ + iy7(p), i(p) = ~ +*y7(p), PeR d , (5.4a) 

where j(p) is a piecewise continuous function of p G M. d with values in and such that 

7(p)p=0, peR d . (5.46) 

One can see that formula (5.3) gives a reconstruction method for Problem 1.1, d > 2, in 
the Born approximation. 

An approximate reconstruction method based on (5.1), (5.2) for Problem 1.1 in di- 
mension d > 2 was proposed for the first time in [Nol] . 
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In the next section we show that, in the Born approximation, Theorem 2.2 (of Section 
2) follows, actually, from (5.3). 

6. Proof of Theorem 2.2 in the Born approximation 

We have that 



vi(x)-v 2 (x)=^J + J je- tpx {vi{p)-V2ip))dp, (6.1) 
\p\<p \p\>p 

|vi(x) - v 2 {x)\ < h{p) + I 2 (p), (6.2a) 



h(p)= J \vi(p)-V2(p)\dp, (6.26) 
\p\<p 

h{p)= J \vi{p) ~ v 2 {p)\dp, (6.2c) 
\p\>p 

where x G R d , p > 0. 

The assumptions ||i>j|| m ,i < R, j = 1, 2, imply that 

Mp)-H P )\< C { f^ , peR d . (6.3) 

Using (5.3) we obtain that 

\vi(p) ~ v 2 (p)\ < C 5 (^)e L 1|^i - $ 2 ||, \p\ < p, 

C 5 (D) = (— ) d [ dx, L = max\x\, ( 6 - 4 ) 

3D 

where || • || is defined according to (2.4). 
Formulas (6.2b), (6.4) imply that 

hip) < C 6 (D)p d e L q^ 1 - $ 2 || < C 6 (D)e Li q^ 1 - $ 2 ||, p > 0, 

C 6 (D) = C 5 (D) J d9, L 1 = L + d. ( 6 - 5 ) 

Formulas (6.2c), (6.3) imply that 

h{p) < C 7 (d,m)Rp-( m - d \ p>0. (6.6) 

Let 

a e]0, 1[ be fixed, 

1-oj (6.7) 
||^>i — ^ 2 || = 5, p = Aln(l+r 1 ), A= — — . 
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Then, due to (6.5), (6.6), 

Hp) < c 6 s(i + (T 1 )^ 1 = c 6 (i + 8y- a 8 a , 
i2(p)<c 7 R(\Hi + s- 1 ))-^ m - d l 

Estimate (2.8) for 6 = ||$i - $ 2 || < 1/2 follows from (6.2a), (6.8). Estimate (2.8) in the 
general case (with modified C2) follows from (2.8) for 6 < 1/2 and the assumptions that 

Thus, in the Born approximation Theorem 2.2 is proved. (This proof is valid also for 
d = 2). 
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